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Abstract 

Motivated by a path planning problem we consider the following procedure. Assume 
that we have two points s and t in the plane and take /C = 0. At each step we add to /C a 
compact convex set that does not contain s nor t. The procedure terminates when the sets 
in K. separate s and t. We show how to add one set to JC in 0(1 + ka(n)) amortized time 
plus the time needed to hnd all sets of V intersecting the newly added set, where n is the 
cardinality oi 1C, fc is the number of sets in K. intersecting the newly added set, and a(-) is 
the inverse of the Ackermann function. 


1 Introduction 

Consider the path planning problem from robotics, also known as the piano mover’s prob¬ 
lem m m Ch.l3]: Given an initial and a target configuration of a robot, the task is to decide 
whether the robot can move from the initial to the target configuration without colliding with 
itself or a surrounding object (and to find such a transformation if it exists). The problem is 
typically tackled by setting up a configuration space X where every robot position is encoded 
as a single point. Then X is partitioned into a free space F C X of allowed configurations and 
its complement F = X \ F denoting configurations that collide with obstacles. The initial and 
final state are denoted by two points s and t in F, and the task is to decide whether s and t 
are in the same path-connected component of F. 

The following approach to solve the path planning problem is discussed by Wang, Chiang 
and Yap HU. Assume for simplicity that the configuration space X is a unit cube in For 
any given subcube, which we call box from now, we can decide whether the box is entirely 
contained in F, entirely contained in F, or both contains points of F and F. We color a box 
green, red, or yellow, respectively, depending on the predicates outcome. Now, starting with 
the entire X, we build a quadtree structure and keep subdividing yellow boxes into 2^^ boxes of 
equal size until one of the following events occur: 

(1) Points s and t lie in green boxes and are connected by a path that lies entirely in green 
boxes. Such a path is a solution to the path planning problem. See Figure left, for an 
illustration. 

(2) Each path from s to t intersects some red square. In this case, no collision-free path from 
s to t can exist, and we say that the red boxes separate s and t. See Figure [T| right, for 
an illustration. 
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The described subdivision strategy is also used for the task of segmentation of digital images; 
see [I] and references therein. In that situation, the approach would decide whether the pixels 
s and t belong to the same connected component of the image. 

How quickly can we decide whether one of the two conditions is satished? Condition (1) can 
be easily checked by union-find m- just create a new element for each new green box and make 
unions to keep together adjacent green boxes, always checking whether the boxes containing s 
and t fall into the same set. That means that the amortized complexity of checking condition 
(1) is almost linear in the number of green boxes produced. For condition (2), the case seems 
less clear - an alternative way of phrasing the condition is to check whether the union of green 
and yellow boxes contains s and t in the same connected component. The union-find approach 
cannot directly be applied because yellow regions might turn into red and, therefore, the area 
covered by the boxes may shrink. In this paper, we discuss how to test the second condition 
in the planar case (d = 2). 



Figure 1; Left: Conhguration space with two (convex) holes. When subdividing the marked 
yellow box according to the dashed lines, s and t become connected. Right: Configuration 
space with an annulus-shaped obstacle. When subdividing the marked yellow box, the union 
of red boxes separates s and t, so no path can exist. 

We consider the following generalization of the problem. We have two points s and t in 
the plane. We get a set K, of compact, convex sets in the plane iteratively, adding the sets one 
by one. Each of the sets added to /C is disjoint from s and t. In the motivating problem, the 
red boxes would be the elements of /C. At the end of the insertion of a new compact convex 
set into /C, we want to know whether JC separates s and t. That is, we want to know whether 
each path from s to t has to intersect some element of 1C. Thus, we want a semi-dynamic data 
structure to store K, that allows the insertion of new elements to /C and decides whether 1C 
separates s and t. 

We show that we can maintain 1C under insertions using a slightly more sophisticated union- 
find approach. The time to insert a new set into K, is the time we need to find all the k 
elements of K. intersecting K^, plus 0{k) union-hnd operations. The idea is based on a classical 
parity argument saying that s and t are separated if and only if we can find a closed curve 
contained in the union of the elements of 1C that is crossed an odd number of times by the line 
segment (. from s to t. We maintain a union-hnd data structure for the sets of 1C and augment 
it by storing additional information about the parity of crossings with the line segment 
Using this additional knowledge, we can quickly decide whether adding a new set to 1C forms 
a cycle that separates s and t, and the information can be maintained under unions and path 
compressions without asymptotic overhead. 

If in the motivating subdivision procedure we always subdivide a largest yellow box, we 
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obtain 0 ( 1 ) time per yellow box and 0 (a(n)) amortized time per red box, where n is the 
number of red boxes and a(-) is the inverse of the Ackermann function. Thus, we obtain the 
same asymptotic behavior for testing conditions ( 1 ) and ( 2 ). 

Roadmap. In Section we discuss a criterion to decide when /C separates s and t in the 
static case. In Section we extend this to the semi-dynamic case, where sets get added to fC. 
In Section]^ we discuss the application to the motivating subdivision procedure. 

Our aim is to provide a self-contained exposition. Some of the arguments are an adaptation 
of Cabello and Giannopoulos [3] to this simpler setting, others can be shorten substantially 
using machinery from Algebraic Topology. 

2 Static connectivity 

Let tC denote a finite family of compact convex sets in the plane, and let K denote their union. 
We use the notation IK = \ K. Let s and t be points in IK. 

The set IK separates s and t if they are in different path-connected components of IK. 
Equivalently, IK separates s and t if each path in the plane from s to t intersects IK. We also 
say that K. separates s and t. 

In the next subsection we discuss a criterion to decide when IK separates s and t. The 
criterion is based on considering all polygonal paths contained in IK, and thus is computationally 
unfeasible. In Subsection 12.21 we discuss how this criterion can be checked in the intersection 
graph of 1C, and thus obtain a discrete version suitable for computations. 

We will consistently use Greek letters vr, 7 , r,... only for (polygonal) curves. 

2.1 Topological criterion for separation 

A polygonal curve tt is generic (with respect to s and t) if vr does not contain s nor t and 
the line segment from s to t does neither contain an endpoint of vr nor a self-intersection 
of TT. We will assume in our discussion that all the polygonal curves are generic. We can 
enforce this assumption making a rotation, so that i is horizontal, and replacing the point s 
by s' = s -|- (0, e), for an infinitesimal e > 0. We always use the same perturbed point s'. Since 
tC is finite, separation of s and t with IK is equivalent to separation of s' and t with IK. The 
computations can then be made using simulation of simplicity [H]. 

We hx I as the line segment joining s' and t. The crossing number of i with a polygonal 
curve TT is the number of intersections of i and vr. We denote by cr 2 {i,'n') the modulo 2 value 
of the crossing number of i and vr. Thus, cr 2 (^, vr) = 1 if and only if the crossing number is 
odd. For the whole paper, any arithmetic involving cr 2 (-, •) is done modulo 2. 

It is important to use always the same perturbed point s'. Then, if a polygonal curve vr 
is the concatenation of tt' and vr", we have cr 2 (I', vr) = cr 2 (£, vr') -|- cr 2 (£, vr"). If we would use 
different perturbed points and the common endpoint of vr' and vr" lies in the line segment st, 
then the inequality does not necessarily hold. 

A polygonal curve vr is closed if its endpoints coincide. It is simple if it does not have 
any self-intersections, except for the common endpoint in the case of closed polygonal paths. 
Note that in the following lemma we do not require simple curves. 

Lemma 1. The set IK separates s and t if and only if there exists a closed polygonal curve vr 
contained in IK such that cz 2 {i, vr) = I. 

Proof. We use the following classical argument, which sometimes is an intermediary step to¬ 
wards a proof of the Jordan’s curve theorem: A simple closed polygonal curve vr separates 
s' and t if and only if i and vr have an odd crossing number. See, for example, Mohar and 
Thomassen |9l Section 2.1] for a formal proof. 
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Assume that K contains a closed polygonal curve vr such that i and vr have an odd cross¬ 
ing number. The curve vr may have self-intersections. If vr is not simple, we can split it at 
self-intersections arbitrarily to obtain simple, closed polygonal paths vri,..., vr^ that have, all 
together, the same image as vr. Since we have 1 = cr 2 (.^, vr) = cr 2 (£, vrj), at least one of the 

curves vTj has 1 = cr 2 (£, vTj). Such a curve vTj separates the endpoints of and thus separates 
s and t. It follows that there is no path in \ vr^ from s to t. Since IK C \ vr*, there is no 
path in IK from s to t. 

Assume that there is no path in IK from s to t. Consider the path-connected component A of 
IK that contains s. Since t is in a different cell of IK and /C is a finite collection of compact, convex 
bodies, there exists a simple closed curve vr contained in the boundary of A that separates s 
and t. We can make shortcuts in vr to obtain a simple closed polygonal curve vr' contained in 
IK that separates s and t. (This can be shown formally using the convexity of the elements of 
JC and the compactness of IK.) The resulting simple polygonal path vr' separates s and t, and 
thus the crossing number of i and vr' is odd. □ 

Lemma 2. Let and be two compact convex sets of JC. For any two generic polygonal 
curves it and vr' contained in U Ky with the same endpoints, we have cT 2 {i, vr) = cr 2 (^, vr'). 

Proof. First note that Ky U Ky does not separate s and t. This can be seen as follows. Let 
be the set of directions. Consider the set of directions of the vectors si, for all x G Ky. Since 
Ky is convex and s ^ Ky, this directions cover less than half of §^. A similar statement holds 
for Ky. It follows that there exists some ray from s to infinity in \ {Ky U Ky). Similarly, 
there exists a ray from t to infinity in \ {Ky U Ky). Those two rays and an extra path far 
enough can be combined to obtain a path from s to f in \ {Ky U Ky). Thus, Ky U Ky does 
not separate s and t. 

Since KyUKy does not separate s and t, Lemma [^implies that any closed path 7 contained 
in Ky U Ky has cr 2 (f', 7 ) = 0. The concatenation of vr and the reverse of vr' is a closed path 
contained in Ky U Ky and therefore cr 2 (f', vr) -|- cr 2 (f', vr') = 0. □ 

2.2 Criterion on the Intersection Graph 

Consider the intersection graph of 1C and denote it by G. Each element Ky G /C is a node of 
G] we will denote the node by v to match standard graph theory notation. There is an edge uv 
in G if and only if Ky and Ky intersect. The graph G is an abstract graph. Next we provide 
a geometric representation. 

For each node u of G choose a point py in Ky. For each edge uv of G, let '){uv) be a 
polygonal path from py to py contained in the union Ky U Ky. Since Ky and Ky are convex 
and intersect, we can always choose ^{uv) with at most 2 edges. The pair 

{{py I V G V{G)},{'y{uv) I uv G E{G)}) 

is a drawing of G. (It is not necessarily an embedding because drawings of edges may cross, 
for example when four axis-parallel squares have disjoint interiors but share a vertex.) For 
each walk W = ei... in G, let '){W) be the polygonal path obtained by concatenating 
7 (ei),..., 7 (efc). If IF is a closed walk, then 7 (IF) is a closed polygonal curve. 

Lemma 3. The set IK separates s and t if and only if there exists a closed wallc W in G such 
that cr 2 (f', 7 (IF)) = 1. 

Proof. Assume that IK separates s and t. Because of Lemma there is some polygonal curve 
vr contained in IK such that cr 2 {i, vr) = 1. We break the path vr into pieces such that each piece 
is contained in the union of 2 sets from JC. Let vri,..., vr/j be the resulting pieces, each of them 
a polygonal curve. For each piece vr*, let Xi and yi be the endpoints of vrj, and let Ky^ and Ky^ 
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be the elements of JC that contain Xi and yt, respectively, so that vr, is contained in Ku. U K^.. 
Note that UiVi is an edge of G. Let W be the closed walk with edges uivi, ..., UkVk- 

We claim that cr 2 (^, 7 (VL)) = cr 2 (t',vr) = 1. To see this, consider for each piece TTj the 
polygonal curve 7 * from to pvi obtained by concatenating the line segment from to Xi, 
followed by vTj, and followed by the line segment from yi to . See Figurej^for an example. For 
each piece TTj, the polygonal curves 7 , and 'y{uiVi) have the same endpoints and are contained in 
the union UK^.. Because of Lemmaj^ we have cr 2 (^, 7 i) = cr 2 (£, 7 (uiUi)). It follows that, 
if we define 7 as the concatenation of 71 ,... , 7 ^, we have cr 2 (£, 7 (IT)) = cr 2 (t', 7 ). Moreover, 
cr 2 (.^, 7 ) = cr 2 (^, tt) because 7 is essentially vr with some spokes connecting Xi to p^., where 
the number of crossings evens out. We conclude that cr 2 (I', 7 (IT)) = cr 2 (^, 7 ) = cr 2 (.^, vr) = 1. 
This finishes one direction of the proof. 



For the other direction, assume that G has a closed walk W such that the crossing number 
of i and 7 (bF) is odd. Since the closed polygonal path 7 (VF) is contained in K by construction, 
Lemma implies that K separates s and t. This proves the other direction. □ 

We extend Lemma to a necessary and sufficient condition for s and t being disconnected 
that involves only a few cycles of G. Let T be any maximal spanning forest of G, that is, T 
contains a spanning tree of each connected component of G. For each edge e of G — E{T), let 
cycle{T,e) be the unique cycle in T + e, and let t(T, e) be the curve 'y{cycle(T, e)). That is, 
T{T,e) is the polygonal curve describing cycle{T,e) in the drawing. 

Lemma 4. Let T be a maximal spanning forest of G. The set K separates s and t if and only 
if there exists some edge e G E{G) \ E(T) such that cr 2 (^, r(T, e)) = 1. 

Proof. The essential idea is to use the so-called cycle space of a graph and the fact that 
{cycle{T, e) | e G E{G) \ E{T)} is a basis. We next provide the details using no background. 

Since we can treat each connected component of G (and thus K) independently, we will 
just assume that G has one connected component. This means that T is a spanning tree of G. 

Fix any node r G V (G) and take the point pr G Kr as a basepoint. For each node v G V (T), 
let T[r, v\ be the simple walk in T from r to u. For each edge uv of G we dehne a closed polygonal 
curve X{uv) as the concatenation of 7 (T[r, u]), 'y{uv), and the reverse of 'y(T[r,v]). Note that 
X{uv) is a closed polygonal path through p^. 

When uv ^ E{T), X{uv) is t{T,uv) with a spoke following 'y{T[r,w]), where w is the last 
common node of T[r,u] and T[r,v]. This implies that 

\/uv G E{G) \ E{T) : cr 2 {£,T{T,uv)) = cr 2 {i, X{uv)). (1) 

When uv G E(T), then X{uv) walks twice the same polygonal curve, and therefore 

Muv G E{T) ■. cr 2 {i,X{uv)) = Q. (2) 
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Assume that the points s and t lie in different path-components of K. Because of Lemma 
there exists some closed walk W inG with cr 2 (^, 7(kL)) = 1- Let uivi, ..., UkVk be the sequence 
of edges in VL, where ui = Vk- Using arithmetic modulo 2 we have 

1 = cr2(^,7(VU)) 

k 

= ^cr2(£,7(ujUi)) 

i=l 
k 

= ^ {cr2{£, 7(r[r, Ui\) + cr2{£, -f{uiVi) + cr2{i, 7(7"[r, Uj])^ 
i=l 
k 

= cr2(£, \{uiVi)), 

i=l 

where in the third equality we have used that each node of W is the endpoint of two consecutive 
edges of VU, which implies that the new terms cancel out. This means that, for some edge UiVi 
of VU, we have cr 2 (t', A(ujUi)) = 1. This edge UiVi cannot be in T because of Q. Therefore 
we have some edge UiVi in E{W), where UiVi ^ E{T), with cr 2 (f, A(ujUj)) = 1. Because of 0 
we have cr 2 {£, T{T,UiVi)) = cr 2 {£, X{uiVi)) = 1. This finishes the proof of one direction of the 
statement. 

For the other direction, assume that there exists some edge e G E{G) \ E{T) such that 
cr 2 (£, t(T, e)) = 1. Taking W = cycle{T,e) and using that r(T, e) = 7 (kF) by definition, this 
means that IT is a closed walk in G with cr 2 (^, 7 (kF)) = 1. It follows from Lemmathat K 
separates s and t. □ 

3 Semi-dynamic connectivity 

In this section we discuss the separation of s and t under the addition of new sets to 1C. We first 
describe a standard union-find data structure because we will build on it. Then we describe the 
setting and the notation we will use. It follows a description of the extension of the union-find 
data structure for our setting. Finally, we describe the data structure, its maintenance, and its 
correctness. 

3.1 Preliminaries: Union-find 

Here we review a standard union-find data structure and some of its properties. See [H Chapter 
21], 0 Chapter 5] or [7] for a comprehensive exposition. 

A union-find data structure represents a disjoint set system supporting the operations 
MakeSet (create a new disjoint set with a single element). Union (merge two sets), and 
Find (return a representative of a given set). We can test whether two elements belong two 
the same set by testing whether the output of Find for those two elements is the same. A 
common realization is to represent each disjoint set by a rooted tree in which each node holds 
one element of the set. The root of the tree holds the representative of the set. Each node has 
a pointer to its parent, while the root points to itself. Then Find simply follows the parent 
pointer until it finds the root of the tree. The union operation merges two trees by making the 
root of one subtree a child of the root of the other. Thus, given two elements, we first locate 
the roots of their corresponding trees calling Find, and then we proceed with the union. 

Two optimizations are commonly used to obtain an efficient realization. Union-by-rank 
determines which root gets merged in a union operation: each root has a rank associated to it, 
in an union we simply make the root of lower rank a child of the root with larger rank, and we 
increase the rank of the root if both roots had the same rank. Path compression makes all nodes 
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found on a search path from a node to its root direct children of the root. For later reference 
and modification, we include pseudocode in Figure Combining these two optimizations, 
each operation has an amortized time complexity of a{n), where n is the number of elements 
in the set system and a(-) is the extremely slow growing inverse Ackermann function. See 
references [H Chapter 21], [7] or m for an analysis of the time complexity. 


Algorithm Find(u) 

1. if u 7 ^ parent(u) then 

2. parent{u) •(— Find {parent (u)) 

3. return parent{u) 


Algorithm Union(u, u) 

1. u ^ Find(u) 

2. V ^ Find(u) 

3. if rank{u) > rank{v) then 

4. parent{v) •(— u 

5. else (* rank{u) < rank{v) *) 

6. parent{u) ■(— v 

7. if rank{u) = rank{v) then 

8. rank{v) •(— rank{v) + 1 


Figure 3: The main two operations in the union-find data structure, u and v are nodes of the 
tree. 


3.2 Setting 

Let s and t be two points in the plane. We have a finite family of convex sets /C, all of them 
disjoint from s and t. Following the previous notation, we denote by K the union of the sets 
in 1C, and by G the intersection graph of /C. 

Consider the addition of a new compact convex set Ku to 1C. We use Anew for the resulting 
set. Knew for the union of its sets, and Gnew for the intersection graph of Anew- 

The analysis of our data structure is based on a maximal spanning forest of the intersection 
graph of the convex sets. The definition of such spanning forest is iterative, as follows. Let 
uvi ,..., uvk be an enumeration of the edges incident to u in Gnew- That is, ..., Ky^, are 
the sets of 1C intersecting the new set Ky. We consider adding the edges uvi,... ,uvk to G 
one by one. We thus define Go as the union of G and a new vertex u for Ky. For each index 
^ < j < k, we define the graph Gj = Gj_i -|- uvj. Note that Gnew = G^. The intermediate 
graphs Gi,..., Gk-i are not intersection graphs of 1C or Anew, but something in between. 

If at the time of adding uvj the vertices u and Vj are already connected in the graph Gj_i, 
then we call uvj a cycle edge. Otherwise, uvj merges two connected components of Gj_i and 
we call it a merge edge. Note that whether an edge is a cycle edge or a merge edge depends 
on the order used in the addition of edges. 

Let T be the maximal spanning forest of G. We define Tq as the union of T and a new 
vertex u for Ky. For each index 1 < j < k we define 


Tj = 


Tj_i if uvj is a cycle edge, 

Tj-i + uvj if uvj is a merge edge. 


It is easy to see by induction that, for each index 1 < j < k, Tj is a maximal spanning forest 
of Gj. We define Tnew as T^. Thus Tnew is a maximal spanning forest of Gnew = G^. 

As it was done in Section \2.2\ for each Ku we choose a point pu in Ku and for each edge 
uv we choose a polygonal curve 'y{uv). These choices are made in the first appearance of the 
node or edge, and remain invariant from there onwards. 
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3.3 Augmented union-find 

We maintain a union-find data structure for the connected components of the graphs Gj. Recall 
that Tj is a maximal spanning forest of Gj. For each node v of Gj, we store a parity bit, 
denoted as parity(v), with the following property: 

• If u is the root of a union-find tree, then parity{v) = 0. 

• If u has parent u) in a union-find tree, then parity{v) = cr 2 {i,Tj[w,v]). That is, we look 
at the parity of the crossing number of i with the polygonal curve from py to py, defined 
by the drawing of Tj. 

For the rest of the paper, any arithmetic involving parity bits is done modulo 2. 

We next argue that the correct parity bits can be maintained in the same complexity as the 
union-find operations, assuming that only certain unions are made. That is clear for MakeSet 
by giving the new node parity 0. 

Consider the Find operation, which changes parent pointers due to path compression. Note 
that the graphs Gj and Tj do not change, but the union-find data structure does. Let u, v, w 
be nodes such that, in the union-find data structure, w is parent of v and v is parent of u. 
Note that 

cT 2 {^,-f{Tj[u,w\)) = cr 2 {(;'^{Tj[u,v])) + c-r 2 {i,'^{Tj[v,w])) = parity (u) + parity (v). 

Therefore, when we update parent{u) •(— w, we just have to set parity (u) •(— parity {u)+parity (v) 
to restore parity(u) to its correct value. 

We can now easily realize the augmented path compression. We define an extended function 
FindExt(m) that, for all nodes v from u to the root r of the tree containing u, sets parent{u) = r 
and updates the value parity(v) accordingly. Pseudocode is given in Figure]^ It easily follows 
by induction that FindExt correctly maintains the parity bit of all elements. 


Algorithm EiNDExT(n) 

1. if n 7 ^ parent(u) then 

2. r Fwb{ parent{u)) 

3. parity{u) •(— parity{u) + parity {parent {u)) 

4. parent{u) ^ r 

5. return parent{u) 


Figure 4: Extended find operation for an element u. 

Finally, we discuss the extension Union to UnionExt. Its arguments are two nodes u 
and Vj such that uvj is a merge edge and the union-find data structure stores the connectivity 
of Gj-i. Since uvj is a merge edge, we have Tj = Tj-i + uvj. This means that the sets 
Ky and Ky. intersect but u and Vj were in different connected components of Gj-i. Like 
before, we first find the roots u and v of their trees using EindExt(-). After this it holds that 
parity{u) = cr 2 {^,^{Tj-i\u,u])), and similarly parity{v) = cr 2 {(-,^{Tj-i\v,v])). 

The walk Tj[u,v\ can be split into Tj-i[u,u], uv, and Tj-i[v,v]. Thus we have 

cr 2 {l,-i{Tj[u,v])) = cT 2 {l,-i{Tj_i[u,u])) + ciC 2 {l,^{uv)) + CT 2 {l,-i{Tj_i[v,v])) 

= parity{u) + cr 2 {l, ^{uv)) + parity{v). 

The last values are either available through parity{-) or computable in constant time. If, for 
example, u gets v as its parent, then we have parity{u) = cT 2 {i, j{Tj[u, v])). The other case is 
similar. We provide the resulting pseudocode in Figure 
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Algorithm UnionExt(m, ?;) 

1. it ^ FlNDEXT(tt) 

2. u ■(— EindExt(u) 

3. b ■>r- parity{u) + parity{v) + cr 2 {£, 'y{uv)) 

4. if rank{u) > rank{v) then 

5. parent{v) ■(— u 

6. parity{v) •(— b 

7. else (* rank{u) < rank{v) *) 

8. parent{u) •(— v 

9. parity (u) ■(— b 

10. if rank{u) = rank{v) then 

11. rank{v) = rank{v) + 1 


Figure 5: Extended union for two nodes u and v. 


The properties of union-find imply that each of the extended operations, UnionExt and 
FindExt, has an amortized complexity of a(n), where n is the cardinality of 1 C. 


3.4 Semi-dynamic data structure 

We now describe the data structure to maintain 1 C. The data structure supports one operation: 
add a new compact convex set Ku to JC and then report whether K, U {Ky} separates s and t. 
We use the notation from Sections 13.21 and 13.31 
The data structure has the following elements: 


an augmented union-find data structure as described in Section 3.3 


• for each element Ky of /C, we store the point py] 


• a semi-dynamic data structure DS{JC) that can find, for the new Ky, all the objects of 
K. that intersect Ky. 


The intersection graph G and the maximal spanning forest T are not kept. They are used only 
for the analysis. 

We next describe how to insert Ky. We use the data structure DS{1C) to find the sets 
Ky^, ..., Ky^_ of JC that intersect Ky. We then insert Ky in the data structure DS{1C) to obtain 
DS{ICyew)- We choose a point py in Ky and create a new node u in the extended union-find 
data structure. 

We then iterate over the edges uui,... ,uvk. We first decide whether the considered edge 
uvj is a merge edge or a cycle edge by checking whether FiNDEXT(rt) and FindExt(uj) return 
the same representative. If uvj is a merge edge, we just call UnionExt( u, Uj) and continue 
with the next step of the filtration. 

Otherwise, uvj is a cycle edge, and we proceed as follows. We want to check whether 
cr2{i,T{Tj,uvj)) = cr2{i,T{Tj-i,uvj)) is 1 or 0. For this, we use that u and v have already 
the same parent because of the calls FiNDEXT(tt) and FindExt(u). If we denote such a 
common parent by r, then 


cr2{i,T{Tj,uvj)) = cr2{£,'y{Tj_i[u,Vj])) + cr2{i,'y{uvj)) 

= cr2{i,jiTj_i[u,r])) + cr2ii,'y{Tj_i[vj,r])) + cr2{£,'y (uvj)) 
= parity{u) + parity{vj) -|- cr2{i,^{uvj)). 
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If ci: 2 {i,T{Tj,uvj)) = 1, then we conclude that Knew separates s and t and we finish the 
algorithm. If cr 2 {i, T{Tj,uvj)) = 0, we proceed to the next edge uvj+i. Pseudocode for the 
insertion of Ku is given in Figure This finishes the description of the algorithm. 


Algorithm Adding to 1C 

1. MAKESET(rt) 

2. parity {u) ■<— 0 

3. for G 1C intersecting do 

4. u ^ FiNDExT(n) 

5. V ^ FindExt(u) 

6. if h 7 ^ h then 

7. UnionExt(u, u) 

8. else {* uv a. cycle edge *) 

9. 6 ■<— parity{u) + parity{v) + cr 2 {i,'y{uv)) 

10. if 6 = 1 then 

11. return “/C U {Ku} separates s and t!!’’ 


Figure 6: Procedure for the addition of Ku- 


It follows from the invariants of the extended union-find discussed in Section 3.3 that we 


are correctly computing the value cr 2 {£, T{Tj,uvj)). If cr 2 {i, T{Tj,uvj)) = 1, then Lemma 
implies that Knew separates s and t. From that point on, we only need to remember that s 
and t are separated. 

If cr 2 ('^, T{Tj,uvj)) = 0, then cr 2 {i, t(T, uvj)) will remain 0 for all future maximal spanning 
forests T. This is so because the maximal spanning forest we maintain is monotone increas¬ 
ing; we only add vertices and edges, but never remove anything. Thus, we never need to 
check cr 2 {i,T(T,uvj)) again later. In particular, if 1C did not separate s and t and we have 
cr 2 {i,T{Tj,uVj)) = 0 for all j, then 


VUU' G E(Gnew) \ ^(Tnew) 


cr2( 


?,r(Tnew, W)) = 0. 


Since Tnew is a maximal spanning forest of Gnew, Lemma [^implies that Knew does not separate 
s and t. 

For each edge uvj we make 2 calls to FindExt, at most one call to UnionExt, and 
additional 0(1) work. This means that for each edge we spend 0(a(n)) amortized time, where 
n is the cardinality of 1C. We also need the time needed to find the elements of 1C intersecting 
the new element Ku- We conclude. 


Theorem 5. Let s and t be two points in the plane. There is a semi-dynamic data structure to 
maintain a family 1C of n compact convex sets in the plane under insertions to decide whether 
1C separates s ean t. The insertion of a new set Ku in K. that intersects k sets of K. takes 
0{l + ka{n)) amortized time, plus the time needed to find the k elements of 1C intersecting Ku- 

Of course, once s and t are separated by 1C, the insertion of each new set can be carried 
out in constant time, since we only need to remember that 1C separates s and t. 


4 Application to dynamic connectivity under subdivision 

We consider now the motivating application discussed in the Introduction for d = 2. 

We have two points s and t inside the unit square X. Initially, the box X is colored yellow. 
In each iteration, we take a largest yellow box, subdivide it into 4 subboxes, and color each of 
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them as red, yellow, or green depending on the outcome of some oracle. The boxes containing s 
or t are always colored yellow or green. We want to know at which point the red boxes separate 
s and t, meaning that each path from s to t contained in the unit square intersects some red 
box. 

Boxes are assumed to contain their boundary, so that any two boxes intersect if their 
boundaries intersect, possibly only at a common vertex. 

For our arguments it is convenient to surround X with 8 red boxes of the same size as 
X. This reduces the problem to finding certain curves within the red region. Without those 
additional squares, we should also consider boundary-to-boundary curves. 

We maintain through the algorithm the intersection graph H of the yellow and red boxes. 
This intersection graph H has one node for each box that is yellow or red, and an edge 
between two nodes whenever the corresponding boxes intersect. The graph H is stored using 
an adjacency list representation [U Chapter 22]. The adjacency list of each vertex is stored as 
a doubly linked list. Moreover, for the appearance of a node v in the adjacency list of u, we 
keep a pointer to the appearance of u in the adjacency list of v. With this, we can perform the 
deletion of a node v in time proportional to its degree. 

When we want to subdivide a yellow box Ku represented by a node u, we can locate its 
set of neighbors N = Nh{u) in the graph H, delete u from the graph, subdivide into four 
boxes, create the at most four new nodes representing the yellow and red boxes arising from 
the subdivision of K^, check for intersection each of them against each of the nodes in N, and 
update the graph H accordingly. All this takes time 0(1 + \N\) time. 

If we always subdivide a largest yellow box, there are at most 12 other boxes intersecting 
it. This means that we can update the intersection graph H of yellow and red boxes in 0(1) 
time. For choosing always a largest yellow box, we can use for example a queue for the yellow 
boxes. Thus, we spend 0(1) time per subdivided yellow box and, for each red box, we get 
its neighboring red boxes in 0(1) time. Using Theorem for the red boxes, and a normal 
union-find for the green boxes, as discussed in the Introduction, we obtain the following result. 

Theorem 6. Consider the subdivision procedure described in the Introduction where we always 
subdivide a largest yellow box. We can perform the subdivision until condition (1) or (2) occurs 
in 0{na{n)) time, where n is the number of subdivisions performed. 

Of course we can also perform the first n steps of the subdivision procedure in 0{na{n)) 
time, and correctly report that neither condition (1) nor (2) hold. 

Acknowledgments We thank Chee Yap for posing to us the problem about connectivity 
under subdivisions. 
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